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Abstract 

We recalculated pion virtual Compton scattering in perturbative QCD 
in this paper. Our calculation avoids some deficiencies in existing liter- 
atures, and treats real Compton scattering as a limit case in which the 
mass of the virtual photon equals to zero. Expressions of the hard scat- 
tering amplitudes from 10 independent diagrams are given explicitly in 
the text. By comparing the effects of different distribution amplitudes 
on the physical observables, we studied the self-consistence of pQCD 
calculation of this problem. 



1 Introduction 

Pion virtual Compton scattering (VCS) 7*7r~ — > 77?" via the reaction eir~ — > e^fir~ 
is observed by SELEX Collaboration at Fermi Lab [l|] for the first time. Although in 
the current available kinematical region, the process can not be predicted precisely 
in perturbative QCD (pQCD), it is possible to observe such processes in the pQCD- 
applicable region with the quick development of experimental techniques. Therefore 
it is meaningful to check the pQCD prediction of this process. 
* Mailing address. 



1 




Figure 1: Reaction (a) eit^ — ► e^fir^ can proceed both through (b)&(c) Bethe- 
Heitler process and through (d) virtual Compton scattering. In experiments, the 
two kinds of process can not be separately detected. When the pion momen- 
tum change is large enough in the process, the amplitude M of the process (d) 
^ir^ — > 77r ± can be factorized as (e) M = J dxdyc/)(x)H(x,y)(f)(y), where H can 
be computed perturbatively, please see fig|| 

For this problem, Tamazouzt calculated a very similar process, 77r ± — > / -f*ir ± ] 
Maina and Torrasso [^] calculated it directly and treated the singular integration 
appearing in it carefully; Li and Coriano [|J calculated it by a rather different way. 
One common problem existing in and 0] is: the authors only directly calculated 
5 diagrams contributing to the unintegrated amplitude but gave no prescriptions 
about how to get the other 15 ones, please see fig.|2| and captions there. Ref . did 
not give its expressions for the unintegrated amplitude but claimed consistence with 
0. As was pointed out by the author of 0, the numerical treatments of have 
some defaults. However, in literature f3[, features that deserve further investigations 
still exist after the revision. For example, a nearly jumping change happened to the 
phase of M LR , please see figure 5 of 0. 

Because of these questions, we decide to recalculate this problem in this paper. 
It will be shown that, all the deficiencies in the literatures will not exist in our 
recalculated results. 
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2 Factorization Theorem of pQCD for Pion VCS 



As indicated in the caption of fig.l, a complete physical process eir~ —>■ c-fir' can 
take place through two ways, Bethe-Heitler process and virtual Compton scattering. 
But we will not calculate such a complete process in this paper, please refer to || 
and ||. We will concentrate on the sub-process 7*7r ± — > 77r ± and use M denoting 
the amplitude of it. Besides its inclusion of real Compton scattering as a limit case, 
it is also meaningful in the preparation for the calculation of the complete process 
e7r~ — > e77r _ . 

Factorization theorem states that for an exclusive process || such as pion VCS 
7*7r ± — > •yii ± , the amplitude of it can be written as the following convolution formula, 

M(p; eq -> p'\ e'q') = J dxdy<fr(x, Q)H(xp; e,q^ yp'\ e', q')(f>{y, Q) , (1) 

where p denotes the momentum of the incoming pion, e and q are the polarization 
vector and momentum of the photon respectively, x denotes one of the momentum 
fraction of the valence quarks in the incoming pion, and that of the other one will 
be denoted by x = 1 — x. The primed variables or y are associated with outgoing 
particles. 

The pion distribution amplitude <p(x,Q) in eq.(|l]) absorbs the long-distance dy- 
namics of M and can be derived by non-perturbative methods 0. Appearing of Q 
in it indicates its evolution with the energy scale. In this paper, instead of consid- 
eration of such evolution 0, ||, we will study the following five phenomenological 
models and their effects on the physics predictions |], [10], [11], [12], [13| , 



b as = VSf n x(l - x), 
W = 1.4706\/3/ 7r x(l - x) exp r 



>hs 



x(l — x) 
5^/^(1 - x)(2x - l) 2 , 

/-, / w „ r 0.07062 
8.8763^/^(1 - x){2x - l) 2 exp[ 



x(l — X) 

P 3 = V3Ux(l - x) [0.6016 - 4.659(2x - l) 2 + 15.52(2x - l) A 



(2) 



with the pion decay constant /„- = 93 MeV and the distribution amplitudes nor- 
malized by / dxcj)(x) = \/ r 3f 7T /6, please see fig. 4 for their explicit shapes. From the 
figure, we can see that, relative to <p as and <p cz , the distribution amplitudes fata and 
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<phs suppress the end point region deeply, while function p3 intensifies both the 
near-end-point region and the center region. 

Contrary to <f)(x), the hard amplitude H(xp;eq — > yp';e'q') in eq.(|lD absorbs 
short-distance dynamics of the amplitude, and can be calculated perturbatively on 
the basis of the diagrams in fig|| With leading Fock state of 7r + (tt~ case is similar), 



k + (p)> = E [\ u \( x P))\^ii x P)) - \u\{xp))\d){xp)) 

where i and i denote the color indices, H can be written as 



1,2,3 



(3) 



H(xp, eq — * yp\ e'q' 
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(4) 



with the coupling constants a e , a s and the charge factor or e u eg suppressed for 
the moment. In eq.(f|), pi, P2 and p% denote the momentum transferred through the 
two fermion and one gluon propagators; explicit expressions of (...7 M t a ...)s depend 
on the details of diagrams. Using identity of SU(3) color group, 



(5) 



and some trace making trick, please see |14J , eq.(f|) can be transformed into the 
following form 



H(x, eq -> y, e'q') = ]T H% 

diag 



{diag.) 



-Vxxyy Tr [(l f *--)diag.l 5 p'(lv-)diag.l 5 l 
diag. 



(6) 



By the usual convention, distribution amplitude <f)(x)<j)(y) will absorb the y/xxyy 
factor, so it will not be included in our later expressions for H ee >. 
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Figure 2: Unintegrated (hard) amplitude H can be calculated on the basis 
of these diagrams. Complete H includes other ten diagrams with the photons 
attaching to different quark lines, and contributions of those diagrams to the 
full amplitude M are equal to the above ones except some charge factors. 



3 Unintegrated Amplitude of j*tt 
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In the center-of-momentum frame of outgoing particles 7 and tt^^, please refer to |J, 
we write all the relevant kinematical variables as follows, 

p'" = ^^(1, - sin 0,0,- cos 0), 



q 



1, sin 9, 0, cos 9), 



p"=(p,0,0,-p), <f =(o;,0,0,p), 

£r = 7/|( ' _1 ' _i >°)> e R = -^(O^ - cos6l > -i,sini 



(7) 



e£ = -=(0,l,-i,0), 



V2 
1 
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1,0,0,1). 



Obviously, 9 denotes the scattering angle of the process. According to parity 
invariance and gauge invariance, we only need to calculate three helicity amplitudes 
for the purpose of computing the amplitude of the complete physical process e7r ± — > 
e77r ± , please see [|] and ||. In order to compare with @, we choose to calculate 
H RR , H LR and H +R , while obtain the other five by the following relations, 
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After introducing the abbreviation c = cos §, s 



sinf, S = (to + p) 2 , v = q^qJS 
and v = 1 — v, we can write all the diagrams contributing to H in more economical 
forms. The results are shown in table 1-3, where we have transformed the diagrams 
with two propagators potentially on shell into some equivalent forms in which only 
one of them can go on shell by the following relation, 

f(x,y) f(x,y) (l 1 \ 



(x — a(y) + ie)(x — b(y) + ie) a(y) — b(y) \x — a(y) + ie x — b(y)+ie) ^ ^ 
This is very important in the numerical integration. Considering the fact that each 
diagram of fig.0 has a companion with its photons attached to different quark lines, 
we include a charge factor for each of the diagrams in the table, where 
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For H RR a9 and Hf R 9 '\ we compared our expressions with those of [fj] in the v — > 
limit (in 0, it is Rb — > 0). Except our consideration of 5 additional independent 
diagrams labelled / — > j, all the other terms, labelled a — > e, coincide with [0]. For 
j^(diag.) ^ - n ^ . g because we employ the convention of || for the virtual 

photon polarization vector, which is different from ||, our expression of it does not 
coincide with that of 0. Our convention is very convenient for future calculation of 
the complete process e7r ± — ► e77r ± . 

In the case of v = 0, by adding all the diagrams in table 2 together, we can get 
a rather simple expression for Hlr, 



,,2 1„ 8c" 2 " 2 

l LR 



h t „ = -(- + -? t: " (if 



'3 3 y x(i- y y 

So, in the v — > limit, the amplitude Hlr is a real number, it has no imaginary 
part. About this point, literature JT5J give a general discussion. It should be notified 



that has an error or misprint in giving its expression for Hlr\ v _,q as Hlr — 
Co(ei — e2) 2 {x — y)c~ 2 s 2 . Obviously, such an unintegrated amplitude will give zero 
amplitude M^r in the integration after multiplied by a symmetric factor cj)(x)4>(y), 
please see eq. ([!]). 

From table 1-3, we can see that, one by one, diagrams in the second row of fig.^§ 
do not equal to those of the first row. From the numerical results of later sections, 
we will be able to see that, as a total, the second row diagrams also do not equal to 
those of the first row. So in this problem, the number of independent diagrams is 10 
instead of 5. Of course, the total number is 20 as we indicate in the caption of fig.@. 

4 Analytical Results of Electron VCS and Quali- 
tative Properties of Pion VCS 

From the aspect of experiencing VCS, unpolarized electrons and pions are similar to 
each other, so we can hope cross sections of VCS on the unpolarized electrons and on 
pions have similar v and 9 dependence. Because electron has no internal structure, 
its VCS cross sections can be get analytically. By the same kinematical variables 
as those of pion VCS, neglecting the mass of the electron, we can get the following 
expressions for electron VCS 7*e — > j e /e, 



da RR ~ v\M RR 



2 
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Figure 3: VCS on the unpolarized electron. 



~|[(l-^) 2 + c 4 )] (12) 
da LR ~ v\M LR \ 2 

~ ^V] (13) 
da + R ~ u|M + /j| 2 

~ 2v 2 s 2 (14) 

Fig.|3| illustrated explicit shape of the v and 6 dependence of the cross sections. 
From the figure, we can easily see that, in both the R — > R and L — > R processes, 
large angle scattering cross sections dominate over the little angle ones. While, in 

the H ► R process, the cross sections depend on the scattering angle weakly. As 

we will indicate in the following, in pion VCS, the same properties of the polarized 
cross sections persist. 



5 Numerical Integrations and Results 

We perform numerical integrations and give results for the cross sections as well as 
corresponding phases for different polarized processes in this section. 

With eqs. ([!]), (@), (@) and table 1-3, and using the following relation [|i~6 



^3^2 = l^da^ = -S^ {M? (15) 
dcos6 2 dt 32vr' 1 ' v ' 



S 



and the principle integration formula, 

f(x,y) f f(x,y*) 

dxdy t"~t = / dxdyP- — r — iir dxdyf(x,y)5(x — a) 

x — a{y) + xe J (x — a) Jo 



/■ 

[ dxdy p!M> = /^p jMz/M + /( )]o (16) 

J x — a J x — a a 



we can get reliable numerical results for the cross sections and corresponding phases 
for different polarization processes. One can see that our treatment of the singular 
point appearing in the numerical integration is the subtraction method. As was in- 
dicated by [|l^, [17| in the similar calculation for proton Compton scattering, different 
numerical treatments of the singular points can give consistent results as long as it 
is applied appropriately. 

We get our numerical integrations by VEGAS program [jn|. As did by ||, we 
take a s = 0.3, a e = 1/137.036 in numerical computations. Noting the fact that, H, 
thus M varies with -4, we show the product S f 7 "'^ instead of d ^ e ' {6 a ) in final results. 

b ' r d cos a cos t) 

First, in fig.(|, using distribution amplitude P 3, we compared the cross sections of 
the process 7| / 7r ± — > 7_rvt ± in the following two cases: (i). the unintegrated amplitude 
includes 10 independent diagrams; (ii). the unintegrated amplitude only includes 5 
independent diagrams (the upper part of fig.^1). Obviously, the cross sections from 
the ten diagram contained amplitudes can not be gained by simply multiplying a 
total factor on those from the five diagram contained ones. 

Second, in fig|| we reconstructed the results of |j] for the cross section S* 3 of 
pion VCS using distribution amplitudes (p p3 . In the right part of this figure, we stud- 
ied the scattering angle dependence of the cross section when v — * 0, by letting v = 
-0.8, -0.1, -0.05, -10~ 2 , -10~ 5 , -10~ 8 in stead of v = -1.0,-0.75,-0.50,-0.25, 
as was done in ||. 

In fig.[], we illustrated the v and 9 dependence of the cross section and corre- 
sponding phases for different polarized processes, with distribution amplitude <p as . 

For the H > R process, the cross section decreases with v —> 0, and it equals to zero 

when v = 0. The scattering angle dependence of it is very similar to that of the elec- 
tron VCS. As in the electron VCS, both for L — > R and for R — > R processes, large 
angle scattering cross sections dominate. The difference is, in the little scattering 
angle regions, the cross section decreases as v — » 0, while in the large angle region, 
the trend reverses. When v — > 0, the phase of Mlr — > 180°, and if we redefine the 
domain of the phase angle, it can be set to 0. 

In fig.0, we compared distribution amplitudes from different models and their 
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effects on the physical observables. We must admit that, relative to those given 
by (pas and (p cz , cross sections given by the end-point-region suppressed distribu- 
tion amplitudes 0^ and <p hs suffer some suppressions. We know that, if the very 
end point region of the distribution amplitudes has very important contributions to 
the cross sections of physical processes, pQCD is non-applicable in calculating this 
problem. Now, our results indicate that this is not the case, so our calculation is 
self-consistence. Of course, to reduce the differences from distribution amplitudes, 
careful treatments of the end point region and consideration of the higher order 
corrections are necessary in further study. 

To see this more clearly, we computed another six polarized cross sections for 
virtual Compton scattering process in fig.§. The upper part of fig|5] corresponds to 
the v = —0.05 case while the down part corresponds to the v = —0.8 case. As in 
fig.[7|, we can see that the end-point-region suppressed distribution amplitudes still 
do not give much-suppressed cross sections at most of the scattering angles. 

6 Conclusions 

We recalculated pion VCS in pQCD in this paper, RCS is treated as a limit case 
in our framework. Comparisons with existing literatures and with electron VCS are 
made in the text. Our study of different distribution amplitudes and their effects on 
the cross sections and corresponding phases of the polarized processes indicates that 
the behavior of the distribution amplitudes at the very end point region does not 
have very strong effects on the physical predictions, but careful treatments of the end 
point region of distribution amplitudes are necessary in the further investigations of 
this problem. 
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Figure 4: Left part, the cross section S 3 ^^ when the hard amplitudes 
include ten independent diagrams; Right part, the same quantity when the 
hard amplitudes only include five independent diagrams. 
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Fi gure 5: Left part, results reconstructed for S 3 ^ ^jifg of ||]. Input photon 
virtuality approaches to zero in the way v = —1.0, —0.75, —0.50, —0.25. Right 
part, our results for the scattering angle dependence of the cross section when 
v goes to zero in the way v = —0.8, —0.1, — 10~ 2 , — 10~ 5 , — 10 -8 . 
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Figure 6: Effects of input photon virtuality on the polarized cross section 
and corresponding phases. The distribution amplitude involved is 4> a s- 

As 

v — > 0, the phase of Mlr — ► 180° asymptotically. 
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Figure 7: Distribution amplitudes and their effects on the real Compton 
scattering process. Relative to 4> as and (p cz , functions cfibhi and 4>hs suppress 
the end point region deeply, but the corresponding cross sections only suffer 
little suppressions. 
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Figure 8: Further study of the effects of distribution amplitudes on the 
virtual Compton scattering. 
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